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The magnetic vortex core in a nanodot can be switched by an alternating transversal magnetic 
field. We propose a simple collective coordinate model which describes comprehensive vortex core 
dynamics, including resonant behavior, weakly nonlinear regimes, and reversal dynamics. A chaotic 
dynamics of the vortex polarity is predicted. All analytical results were confirmed by micromagnetic 
simulations. 
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I. INTRODUCTION 

Manipulation of complex magnetization configurations 
at the scales of nanometers and picoseconds is crucial 
for the physics of nanomagnetisrrP. Among the variety 
of different topologically nontrivial configurations special 
interest attracts the vortex configuration: it can form a 
ground state of the micro- and nanosized disk-shaped 
particles (nanodisks). A magnetic vortex is character- 
ized by an in-plane curling flux-closed structure, which 
minimizes the magnetostatic energy of the particle, and 
the out-of-plane region of the vortex core with about 
the size of the exchange length (typically about 10 nm 
for magnetically soft materials-), which appears due to 
the dominant role of the exchange interaction inside the 
cord^. The direction of the vortex core magnetization, 
the so-called vortex polarity p = ±1 (up or down), can be 
considered as a bit of information in the nonvolatile mag- 
netic vortex random-access memories (VRAMj^El. To 
realize the concept of VRAM one needs to control the 
vortex polarity switching process in a fast way. 

There exist different ways to switch the vortex polarity. 
One can distinguish two basic scenarios of the switch- 
ing: (i) Axially-asymmetric switching occurs, e.g. un- 
der the action of different in-plane AC magnetic fields 
or by a spin polarized current, see Ref. [6]and references 
therein. Such a switching occurs due to the nonlinear res- 
onance in the system of certain magnon modes with non- 
linear coupling^, which is accompanied by the tempo- 
rary creation and annihilation of vortex-antivortex pairs, 
(ii) The axially-symmetric (or punch-through) switch- 
ing occurs, e.g. under the influence of a DC transversal 
fielcpHUl^ The mechanism of such a switching is the di- 
rect pumping of axially-symmetric magnon modes. Very 
recently the resonant pumping of such modes by an AC 
transversal field was pr oposed to switch the vortex in mi- 
cromagnetic simulationsISIHl^ which gives a possibility to 
achieve a switching at much lower field intensities. 

The aim of the current study is to develop a theory 
for the axially-symmetric vortex polarity switching. We 
propose a simple analytical two-parameter cutoff model, 
which allows to describe the main features of the compli- 



cated vortex dynamics under the action of AC pumping, 
including nonlinear resonance and magnetization rever- 
sal. Our model predicts a chaotic dynamics of the vortex 
polarity, which is analyzed in terms of Poincare maps. 
Our full-scale micromagnetic simulations confirmed all 
analytical predictions. 



II. TWO-PARAMETER CUTOFF MODEL 

We consider the model of a classical 2D Heisenberg fer- 
romagnet with effective easy-plane anisotropy, caused by 
the dipolar interaction, under the action of a transversal 
AC field. The energy of such a magnet, normalized by 
the value ttA with A being the exchange constant reads: 
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Here m and (j) ^re related to the compo- 
nents of the magnetization vector M = 
Ms (Vl — 'm? cos 0, Vl — Tn? sin 0, m) , the parame- 
ter i = ^ Aj^-nM'g is the exchange length, Ms is 
the saturation magnetization, and hir) — hsinojT is 
the dimensionless external AC field. We use here the 
dimensionless time r = riot with Qq — 4:njQMs with 
7o being the gyromagnetic ratio. The magnetization 
dynamics follows the Landau-Lifshitz equations, which 
can be derived from the following Lagrangian 
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Here and below the overdot means the derivative with 
respect to r, the parameter rj is the Gilbert damping 
coefficient. 
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In order to describe the switching phenomena we pro- 
pose a simple analytical picture using the following two- 
parameter Ansatz for the magnetization variables: 

m(r,r) =/i(T)/Q) , 0(r,T) ^xil+V-MgQ) . (4) 

We consider the vortex core amplitude m(0,r) — [lij)^ 
which direction has the sense of the dynamical vortex 
polarity and is considered as a collective variable to- 
gether with the in-plane turning phase 'ipir), the func- 
tions f{x) and g{x) describe the vortex structure. We 
use a Gaussian distribution for both functions, f{x) — 
g{x) = exp(— a;'^/2), which is in a good agreement with 
simulation data. One has to note that such an Ansatz 
describes an axially-symmetric vortex solution together 
with the simplest axially-symmetric magnon mode: the 
real solution just slightly varies the profile of the func- 
tions f{x) and g{x). Also it is possible to take into ac- 
count higher modes (with additional nodes on r), but we 
try to make the picture as simple as possible. 

Using this Ansatz one can calculate the total energy of 
the vortex state disk as follows: 
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FIG. 1: Effective double-well potential see ^ 

(solid line). The cutoff parameter A = 0.07. The typical 
evolution of the energy as function of the dynamical 

polarity /Lt(r) under the periodic pumping for 
unidirectional switching is shown by the thin line as a 
result of numerical integration of Eqs.(|8]). System 
parametersp: h = 0.001 with frequency uj = 0.7547, 
damping rj = 0.001, initial conditions 
^(r = 0) = /io = 0.384, ■0(t = 0) = 0, integration time 
Tall = 5000, the trajectory is shown till Tgh ~ 500. 



Here R is the disk radius and a is a cutoff parameter, 
which is of order of the magnetic lattice constant ag . We 
introduce here the cutoff in order to take into account 
discreteness effects. It is worth to remind that in con- 
tinuum theory two vortex states with different polarities 
are separated by an infinite barrier which prohibits the 
switching in a simply connected domain. The function 
(o'(/i, ^p) is the effective energy of the system: 
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which finally read: 
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We start with the case without damping, 77 = 0. In 
this case one can easily exclude the turning phase ij) from 
the consideration, which results in the following effective 
Lagrangian for the dynamical polarity only: 



Here Li2 {x) is a dilogarithm function^^, >r = lnA + 7+ l 
with 7 « 0.577 being Euler's constant, and we assume 
that A = aV^^ < 1. 

The effective Lagrangian and dissipation functions take 
the following forms: 
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The effective equations of motion are then obtained as 
Euler-Lagrange equations 
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The Lagrangian ^ describes the motion of a particle 
with the variable mass ^{^) in the double-well potential 
(/i) under the action of periodical pumping. The typi- 
cal shape of the potential (/z) is shown in Fig.[l] it has 
two energetically equivalent ground states with ^ = ±/io, 
which correspond to vortices with opposite polarities. In 
our cutoff model /zq is a nonzero solution of the transcen- 
dent equation: 
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dXrTrWrW.^ ^^ = {M>}, (7) 



For A = 0.07, the energy minimum corresponds to /ip — 
0.384, see Fig. [l] One has to note that our model works 
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FIG. 2: The nonlinear resonance: amplitude-frequency 
characteristics from numerical solution of 
Eqs. ^ (dashed curve) with initial conditions 
fi{T = 0) = and tplr = 0) = 0. The cutoff parameter 
A = 0.07, the damping t] ~ 0.001. The thin curves 



correspond to the analytical solution ( 14 ) without 



damping: the solid lines correspond to the stable 
solution and the dash-and-dot line to the unstable 

The field amplitude h/rj = 0.15. The total 



region (15 



computation time for each point Ttot = 3000. 



FIG. 3: Diagram of dynamical regimes for the field 
parameters (amplitude and frequency) as numerical 
solution of Eqs. ([s]): circles correspond to one-period 
oscillations, triangles correspond to multiple-periods 
oscillations, and squares correspond to chaotic 
dynamics. White domains correspond to the absence of 

switching. Parameters are the same as in Fig. [2] (the 
total computation time Ttot depends on the regime, see 
the main text). 



only for < 1. Another method is to work with a 
trigonometric variable, the vortex core angle instead 
of the dynamical polarity fi = cosi?. We checked that 
the usage of "i? provides the same physical picture, but 
the effective equations look awkward, so we keep to work 
with /i. 

The dynamical polarity ij.{t) satisfies the following 
equation, see ([9]): 



xfj, ln(l - /i^) _ 2h . 
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The linear oscillations near the potential well bottom 
have the usual harmonic shape: 



with = 1/(4 — /ip) being the effective mass of a small 
oscillating particle near the well bottom and ujq being the 
eigenfrequency. 

Let us study the weakly nonlinear dynamics using the 
method of multiple scaleJiStll!. We limit ourselves by the 
three-scale expansion as follows: 



oj = uiQ + e 0J2, /i = e <C 1, 



(13) 



which provides a valid weakly nonlinear expansion under 
condition that the field amplitude is much less than the 



frequency detuning, h/{uj — ujq) <^ 1. Eq. (Ill together 
with an expansion ([13^ results in the set of equations for 
/i„, see Appendix |A| for details. Following the Floquet 
theorjfi^ one has to remove the mixed-secular terms in 
such equations, which finally provides the nonlinear res- 
onant curve (see Appendix [A| for details): 



w±(|a|) = cjo - ci\a\ ± C2t 



(14) 



where \a\ is an amplitude of oscillations, see ( |12[ ), the 
parameters ci and C2 are calculated in ( A3 ) . A typical 



nonlinear resonance curve is plotted in the Fig. [2| The 
low frequency branch contains a shock-stalling region. 
The upper limit of such an instability region can be found 
using the condition dtu^/da — 0, which finally results in 
the limit frequency 
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The unstable part of the resonance curve is plotted 
in Fig. [2] by the dash-and-dot line. Further increase of 
the field amplitude leads to a broader instability domain. 
Moreover, we will see below that stronger pumping re- 
sults in an essentially different kind of dynamics, leading 
to the switching of the vortex polarity and to chaotic 
behaviour. 

If we increase the amplitude of the forcing, the system 
goes to the strongly nonlinear regime. We analyse such 
regimes using numerical solutions of Eqs. ([s]). First of 
all, the regular oscillations of the dynamical polarity /^(r) 
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(a) h/rj = 11, uj/ujQ = 1.4 , 15 000 points 
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FIG. 4: Poincare map for regular oscillations. The circle 
corresponds to h/i] = 5, uj/ujo — 0.7, triangles 
correspond to h/rj ~ 11, cj/wq ~ 2, and squares 
correspond to h/rj — 13, lj/ljq = 1.95. Parameters are 
the same as in Fig. [2j 



between two potential wells occur in a wide range of pa- 
rameters (the typical oscillations (fj,) are plotted by the 
thin curve in Figjn. The dia gram of dynamical regimes 
is shown in Fig. 131 Different types of dynamical regimes 
are classified in accordance to the Poincare maps. These 
maps are constructed for 20 000 periods of the field oscil- 
lations for frequencies greater than O.Qujq and for 15 000 
oscillation periods for other frequencies. The first 5000 
points are dropped from consideration in order to exclude 
transient processes. The diagram of dynamical regimes 
has a resonant behaviour at the frequencies uJo/3, wo/2, 
ujo and 2wo, see Fig.jSj 

There are three different dynamical regimes on the di- 
agram: (i) The one-period oscillations (circles in Fig. [s]) 
occur in a wide range of parameters, generally in the 
vicinity of the resonance frequency wq, the Poincare map 
for this regime has one stable focus (a circle in Fig. l4|. 
(ii) The multiple-period oscillations (triangles in Fig.Tsl) 
occur typically near the doubled resonance frequency; the 
Poincare map has a few points which are attended every 
pumping period and the trajectory in phase space (/x, "0) 
makes a few windings before closing, see Fig. [4] 

(iii) The chaotic oscillations of the dynamical polarity 
/J, (squares in the Fig. [s]) take place in the transition re- 
gion between the oscillations of the types (i) and (ii) . The 
corresponding Poincare map has the shape of a strange 
attractor, see Fig. [5a] Apart them the chaotic dynamics 
occurs at the resonance frequency the weak enough field 
amplitude and in the wide range of the low frequency 
pumping. The low-frequency dynamics also corresponds 
to a strange attractor, a typical picture is presented in 
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FIG. 5: Strange attractors on the Poincare map. 
Equilibrium polarities are marked by the circles. 



III. NUMERICAL STUDY OF THE DIFFERENT 
DYNAMICAL REGIMES 

In order to check all predictions of the two-parameter 
cutoff model, we performed a full-scale numerical mod- 
elling using the OOMMF frameworlP^l^ which simulates 
the Landau-Lifshitz equations. Numerically we modelled 
a cylinder-shaped sample with radius R = 99 nm and 
height L = 21 nm using material parameters for Permal- 
loy (Ni8oFe2o): exchange constant A = 26pJ/m, satura- 
tion magnetization Alg — 860kA/m and Gilbert damp- 
ing coefficient a = 0.01. The two-dimensional space 
mesh 3 X 3 X 21 nm is used. Initially, the vortex has 
an upward polarity and a counter-clockwise chirality. 

The conditions of our numerical experiment were simi- 
lar to simulations by Wang and Dongl^ and Yoo et al^. 
However, our task was to check the new dynamical 
regime of chaotic dynamics. That is why we need to 
study the long time dynamics. 

First we examined the resonant frequency of radial spin 
waves. A 30 mT constant pulse during 100 ps was ap- 
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FIG. 6: The instability domain near the first resonance 
frequency for simulations with the field amplitude 5 mT. 
The region near the first axially-symmetric harmonic is 
shown. 



plied to the sample; it excited low-amplitude spin waves. 
By analysis of a Fourier spectrum for a 3.7 ns long time 
dynamics of the total magnetization along the disk axis 
we calculated the eigenfrequency of the lowest spin wave 
mode equals /o = 14 GHz. 

We study the vortex dynamics under the action of the 
sinusoidal magnetic field 

B{t) = e^B sm2Tr ft, (16) 

directed perpendicular to the face surfaces of the sample. 
The vortex dynamics under such a field has a resonant 
behaviour. The weak pumping causes the resonance on 
the frequency /q. If we increase the field amplitude, the 
system goes to the nonlinear regime. The weakly non- 
linear regime corresponds to the nonlinear resonance, see 
Fig. (the resonance on the first axially-symmetric har- 
monic). 

To systemize the complicated dynamics of the vortex 
polarity, we compute the phase diagram of the switch- 
ing events by varying the field frequency in the vicin- 
ity of /o from 9 to 19 GHz with steps of 1 GHz and the 
field amplitude from 10 to 110 mT, see Fig. [7] There are 
two strong res onanc es in this range, which agree with 
previous resultJ^^Ell. The lower resonance frequency is 
located between 12 and 14 GHz; it corresponds to the 
axially symmetric mode without radial nodes (m = 0, 
n = 0). The second resonance is located near 18 GHz; it 
corresponds to the axially symmetric mode with a single 
radial node (m = 0, n = 1). Since we expect a chaotic 
dynamics, it is necessary to analyse the long-time be- 
haviour: numerically we checked the magnetization state 
every picosecond during 10 ns interval. The chaotic vor- 
tex polarity dynamics during this time is observed in 14 
simulations (see filled squares in the Fig. [7^) where the 
vortex polarity switching mechanism corresponds to the 
axial-symmetric way. The typical shape of oscillations is 



presented in the Fig. [7]d) for i? = 70 mT and / = 14 GHz. 
We examine the character of the polarity oscillations by 
an autocorrelation function 

^(^^) = ^EMi.+jXij), hJ=^, (17) 

3 

where the discretized time t„ with steps At = t„+i — i„ = 
1 ps are used, and N = is a number of snapshots. The 
function fJ.{ti) is the discrete dynamical polarity, which is 
defined as average magnetization of four cells in the cen- 
ter of the vortex core, normalized by the magnetization 
in the absence of the forcing. For the chaotic signal C (t) 
rapidly decays, see Fig. [Tj;). We marked points on the 
diagram of switching events by filled squares for simula- 
tions where the autocorrelation function rapidly decays 
and the distance between the maximum of the autocor- 
relation function and the first zero is smaller than 1 ns. 
Plots of C (t) for other simulations with chaotic dynamics 
look similarly. 

In all simulations the set of first switchings occurs dur- 
ing the first nanosecond and is accompanied by a high- 
amplitude axially-symmetric spin wave radiation. How- 
ever, typically, the vortex position at the origin is unsta- 
ble; during the field pumping the higher axially nonsym- 
metric modes (m ^ 0) can be excited, which causes a 
vortex motion towards the disk edge surface. In such a 
case the switching occurs through the axially- asymmetric 
mechanism, which is accompanied by the temporary cre- 
ation and annihilation of a vortex-antivortex pair, see 
Ref. [B]and references therein. Such switching events are 
shown in the Fig. [7^) by the filled diamonds. We do not 
analyse them due to an insufficiently short time interval, 
compared to the relaxation time, which corresponds to 
the axial-symmetric switching scenario, discussed in this 
work. 



IV. CONCLUSIONS 

The axially-symmetric vortex polarity switching is an 
efficient way for the magnetization reversal on a sub- 
nanosecond time scale. Very recently such a scheme was 
realized by the micromagnetic simulations in Refs.[T51and 
[Ml To gain some insight to the resonant switching effect, 
Wang and Dong computed an exchange field inside the 
vortex core: it changes rapidly during the vortex rever- 
sal. Yoo et al.^ noticed that the switching occurs only 
if the exchange energy exceeds a threshold value. The 
crucial role of the exchange interaction becomes clear in 
the analytical approach developed in the current study. 
Our two-parameter cutoff model explains the switching 
phenomenon in terms of the nonlinear resonance in a 
double-well potential. Such a potential arises mainly 
from the exchange interaction: the presence of two wells 
corresponds to the energy degeneracy with respect to the 
direction of the vortex polarity (up or down) ; the energy 
barrier between the wells becomes higher as the discrete- 
ness effects become less important. 
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a) Diagram of switching events 
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FIG. 7: Numerical simulations of the vortex dynamics: a) Diagram of dynamical regimes for the field parameters 
(amplitude and frequency). Diamonds and squares correspond to the switching events (the filled diamonds mark 
parameters, where the vortex rapidly moves away from the origin after a few switchings, the open diamonds mark 
the vortex which does not change its position) ; the filled squares show parameters for the chaotic polarity dynamics 
and the crosses correspond to the region without switching events, b) Example of the chaotic dynamics for the 
applied field amplitude B — 70 mT and the frequency / — 14 GHz. c) Autocorrelation function for the plot b). 



In terms of our model the switching can be considered 
as the motion of an effective mechanical particle with a 
variable mass in the double-well potential. Under the ac- 
tion of periodical pumping the particle starts to oscillate 
near the bottom of one of the wells. When the pumping 
increases, there appear nonlinear oscillations of the par- 
ticle; under a further forcing the particle overcomes the 
barrier, which corresponds to the magnetization reversal 
process. The chaotic dynamics of the magnetization is 
an analogue of the chaotic oscillations, e.g., in a Duffing 
oscillator^. 

In summary, we analyse analytically and numerically 
the axially-symmetric scenario of the vortex polarity 
switching, induced by an alternating magnetic field di- 
rected perpendicular to the nanodot surface. We propose 
a simple analytical two-parameter cutoff model, which 
describes the vortex polarity dynamics under such a reso- 
nance pumping and shows the possibility of both periodic 
and chaotic polarity oscillations by Poincare maps. The 
micromagnetic simulations for Permalloy confirm a vari- 
ety of the dynamical regimes and confirm our analytical 
predictions. 
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Appendix A: Analysis by the Method of Multiple 

Scales 



We use the method of multiple scalea^SHlSI to treat an- 
alytically Eq. (11 1. We limit ourselves to the three-scale 



expansion (13 1. Since we have three different time scales 



To, Ti, and T2, one has to modify the time derivatives as 
follows: 



n=0 



d 

d7^' 



(Al) 



The equations governing /j,2, and ^3 are 



(A2a) 



(A2b) 



+ -00132^^1 + DoDifi2] } + h3 sm{LJoTo + U2T2), (A2c) 
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where we used the following notations: 



ki 



-9A.g-M^ + ><(l-Mg)'(4 + 3/.g) 

Mo (1 - Mo) 
3(12 - 27^Il + 38A^g + M^) + 3>r (1 - /xg) 



The solution of the Eq. (A2a) reads /ii — 
y4(Ti,r2)e'""^" + A*(ri,T2)e-*^^ To prevent 



the secular terms in the Eq. (A2b), one has to put 
A{Ti,T2) = A{T2); the same condition for Eq. ( |A2c 
gives an equation for the oscillation amplitude of fii. 



D2A{T2)+4iciA^A* = -h3^oC2e"^'^' 
Sfcoifci 3fc2 5fcJ 



Cl 



C2 



1 



(A3) 



By solving the Eq. (A3) one gets finally the Eq. (f4). 
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